Non-thermal hot electrons ultrafastly generating hot optical phonons in graphite 
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Ultrafast dynamics of graphite is investigated by time-resolved photoemission spectroscopy. We 
observe spectral features of direct photoexcitations, non-thermal electron distributions, and recovery 
dynamics occurring with two time scales having distinct pump-power dependences. Additionally, 
we find an anomalous increase of the spectral intensity around the Fermi level, and we attribute 
this to spectral broadenings due to coupled optical phonons in the transient. The fingerprints of the 
coupled optical phonons occur from the temporal region where the electronic temperature is still 
not definable. This implies that there is a mechanism of ultrafast-and-efficient phonon generations 
beyond a two-temperature model. 



Investigation of the non-equilibrium dynam- 
ics after an impulsive impact provides insights 
into couplings among various excitations. A 
two-temperature model (TTM) is often a start- 
ing point to understand the coupled dynamics 
of electrons and lattice vibrations: the optical 
pulse primarily raises the electronic tempera- 
ture T e i while leaving the lattice temperature 2] 
low; subsequently the hot electrons heat up the 
lattice until T e i = TJ is reached. This tempo- 
ral hierarchy owes to the assumption that the 
electron-electron scattering rate is much larger 
than the electron-phonon scattering rate. We 
report herein that the TTM scheme is seri- 
ously invalidated in semimetal graphite. Time- 
resolved photoemission spectroscopy (TrPES) of 
graphite reveals that fingerprints of coupled op- 
tical phonons (COPs) occur from the initial mo- 
ments where T e \ is still not definable. Our study 
shows that ultrafast-and-efficient phonon genera- 
tions occur beyond the TTM scheme, presumably 
associated to the long duration of the non-thermal 
electrons in graphite. 

One of the most interesting observation in the ultrafast 
dynamics of graphite is that T e i in the sub-picosecond 
temporal region stays somewhat low even though irradi- 
ated with an intense femtosecond optical pulse This 
led to the picture that the electronic energy is quasi- 
instantancously transferred to the COPs through strong 
electron-phonon couplings [l| based on the TTM scheme 
However, electron-phonon coupling constant is 
reported to be moderately small in graphite [fjj|, mak- 
ing the mechanism of the ultrafast COP generation elu- 
sive. Subsequent studies also suggest the nearly instanta- 
neous COP gene ration coupled to the electron dynamics 
@, S, !, EE E, E Q, which is considered to af- 
fect ballistic transports at high fields [U, El E, E E] ■ 
Nevertheless, direct observation of the electron distribu- 
tion in the transient is limited 
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FIG. 1: TrPES system, (a) A schematic of the TrPES system, 
(b) Spectra of graphite recorded in normal-emission geometry. 
The increase of the intensity at Ef with increasing T is due to 
the increased population of the COPs [f| . Inset of (b) shows 
the surface Brillouin zone of graphite, and the area probed in 
the normal-emission geometry is indicated by a circle. 



over, simultaneous detection of the electron distribution 
and the phonons in the transient has been beyond reach. 
TrPES [Fig. Ufa)] is one of the most powerful tools to 
investigate the dynamics of the electrons, since it can 
provide information of the transient electron distribu- 
tions in a wide e nerg y range across the Fermi level {Ep) 
0, E3,|2ll Hi| . Furthermore, Liu ef al. 0| recently 

reported that fingerprints of COPs occur in the photoe- 
mission spectra of graphite, as we shall explain below. 
Therefore, it became possible to monitor simultaneously 
the electron distribution and the fingerprints of COPs 
during the ultrafast dynamics of graphite by TrPES. 

The fingerprints of COPs show up in the spectra 
recorded in a normal-emission geometry [f|, that is, when 
we detect the photoelectrons around V of the surface 
Brillouin zone of graphite [inset in Fig. HI))]. Since 
there are no bands around T in the vicinity of Ep, the 
signal consists of photoelectrons around K (K') indi- 
rectly scattered into the vicinity of T mainly by phonons. 
As we shall see later, signals of direct photoexcitations 
around the K (K r ) point indeed occur in the TrPES spec- 
tra recorded in the normal-emission geometry. A gap- 
like feature of ~70meV occurs in the spectra near- Ef 
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FIG. 2: TrPES of graphite Spectra recorded at 17-mW pump at t < ps (a, b), < t < 0.2 ps (c, d), and t > 0.2 ps (e, f). 
Here, (b, d, f) are semi-logarithmic plots of (a, b, c), respectively. Insets in (d) and (f) show spectra at t = ps recorded at 
various pump powers and temperatures, respectively, (g) Spectral weight between 0.8- 1.0 eV as a function of t overlaid with 
a Gaussian, (h) Spectral weight under 0.75-1.8 eV (above the cutoff; circle) and under 0.0-1.8 eV (in the unoccupied side; 
square) at t = ps as a function of pump power overlaid with first and second order polynomial functions, respectively. The 
schematic in (c) shows direct excitations from the n bands to the n* bands. 



at low T [Fig. QIb)], since phonon absorption process 
is quenched and phonon emission (energy-loss) process 
dominates. With increasing T, phonon absorption pro- 
cess is increased, and the spectral weight tails into higher- 
energies resulting in an increase of the spectral intensity 
at Ep [Fig. H£b)]. This is similar to anti-Stokes lines in 
Raman spectra gaining stronger intensity at higher tem- 
peratures. Thus, we utilize the spectral intensity at Ep 
in TrPES as a measure of the number of COPs in the 
transient. The optical phonons monitored herein are as- 
signed to the 67-meV out-of-plane and the ~150-meV 
in-plane modes [f|. 

Results. — Figure Eta-f) show TrPES spectra of 
graphite, I(uj,t), recorded under a pump power p=17 
mW (a fluence of ^14 /iJ/cm~ 2 ) at room temperature. 
A movie file is provided in Supplementary Information. 
Overall, we observe that the electrons are pumped from 
the occupied side to the unoccupied side and subsequent 
recovery dynamics lasting over several tens of picosec- 
onds: at t < ps, the spectral intensity I(uj,t) is in- 
creased in the unoccupied side; at < t < 0.20 ps, I(lu 7 t) 
at ui ~ 0.75 eV is decreased, whereas that at ui ~ eV is 
increased; at t > 0.20 ps, I(uj, t) in the unoccupied side is 
decreased. Here, t = ps has been determined utilizing 
the fast response of ^20 fs observed at lu > 0.8 eV [Hj], 
and the time resolution [full width at half the maximum 
(FWHM)] is estimated to be zii = 0.43ps, see Fig.^g). 



We observe a plateau feature in the unoccupied side 
during the pump [Fig. [2]Jb) and the movie file in Sup- 
plementary Information], which is a hallmark of a non- 
thermal electron distribution [22[ ■ The plateau turns over 
into an exponential tailing at t > 0.20 ps [Fig. [2Jf)], in- 
dicating that the electrons are distributed according to 
the Fermi-Dirac function. That is, electronic thermaliza- 
tion occurs at r e ~ 0.2 ps so that T e ; becomes definable 
thereafter. In metallic materials, typical time scale for 
electronic thermalization is considered to be r e ~ 10 fs or 
less Q , and if At^>r e , one would not expect to observe 
a non-thermal distribution of the hot electrons. In fact, 
in the TrPES study of a metallic Bi 2 Sr 2 CaCu20 8+( 5 0], 
Perfetti et al. observed with At ~ 90 fs that the spectra 
mostly obey Fermi-Dirac statistics even at £~0ps. The 
smallness of l/r e in graphite can be attributed to the 
semimetallic band structure: Since electronic states near 
Ep occur only around K(K') points of the surface Bril- 
louin zone, the available phase space for electron-electron 
scatterings becomes vanishingly small in approaching Ep 
particularly after the initial avalanche of the hot electrons 
towards Ep. This can act as a bottleneck for the non- 
thermal electrons to relax into the thermal (Fermi-Dirac) 
distribution. 

The plateau observed just after the pump extends up 
to a cutoff at w~0.75eV, as indicated by a bar in Fig. 
&h). This cutoff can be understood as a fingerprint of di- 
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rect photoexcitations occurring around the K(K') point: 
since the n- and 7r*-bands are nearly symmetric about 
Ep as shown in the schematic in Fig. [2J direct excita- 
tions are dominated by the transitions from uj = —hvi/2 
to +hvi/2, and therefore, the cutoff energy can be identi- 
fied to u ~ hui/2 = 0.75 eV. Note that such a fingerprint 
of direct excitations would not be detected if r e <C At. 
Therefore the ~0.75-eV cutoff strengthens the conclu- 
sion that the duration of the non-thermal electron distri- 
bution is detected in the present study. 

With increasing p, the cutoff is blurred [inset in Fig. 
[2jd)] , so that the unoccupied side of the spectra becomes 
featureless and the line shape becomes similar to the TY- 
PES spectra of graphite reported previously |21[. The 
intensity above the cutoff at t — ps grows quadrati- 
cally with p while the intensity in the unoccupied side 
grows linearly with p [Fig. GJh)]. The latter indicates 
that the number of electrons excited by the pump is pro- 
portional to p at least up to 78 mW, whereas the former 
indicates that the excited electrons are further scattered 
above the cutoff by other excitations such as the hot elec- 
trons themselves or hot phonons, since such scatterings 
occur roughly proportional to the square of the popu- 
lation of the excited particles. The intensity and the 
sharpness of the cutoff are almost independent of T [in- 
set of Fig.[2Jf)], indicating that the blurring of the cutoff 
is not related to the thermally populated excitations in 
the initial state. 

We now turn to the variation of the spectral weight 
around Ep, which serves as the fingerprints of COPs in 
the transient, as explained previously. One can see that 
the variation starts from the beginning of the transient in 
accord with the reports of the nearly instantaneous gen- 
eration of the COPs in graphite [1, [HI, [3] ■ Our findings 
therefore show that the ultrafast COP generation takes 
place from the temporal region where the hot electrons 
are not thermalized [see, Fig. [Ha) and the movie file in 
Supplementary Information]. After t ~ 0.2 ps, the in- 
tensity at Ep starts to decrease, see Fig. [2{e), indicating 
that the COP generation is mosltly accomplished within 
t < 0.2 ps. The results are in strong contrast to the TTM 
scheme, where COPs are assumed to be cool at the be- 
ginning and then gradually heated up by the electrons 
that are thermalized. 

Analysis. — First, we investigate how the total elec- 
tronic energy dissipates with time. As I(uj,t) reflects 
the occupied density of states (DOS), J °° u)AI(w,t) duj 
is a measure of the excess electronic energy, and we plot 
this as a function of t in Fig. EJa). Here, AI(uj,t) = 
I(uj,t) — I(u,—5ps). The energy dissipation occurs 
with two time scales having distinct pump-power de- 
pendences: at t < 1 ps, the energy dissipation rate pos- 
itively depends on p, whereas at t > 1 ps, it is indepen- 
dent of p. The fast dynamics at t < 1 ps is attributed 
to the net energy flow from the hot electrons to the 
hot COPs through electron-phonon scatterings, since this 
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FIG. 3: Decay rates, (a) Electronic energy (see text) vs t 
for various pump powers. Each curve has an arbitrary offset. 
Dashed line represents the time resolution, (b) The spectra of 
decay time and decay rate (left and right axes, respectively) 
of the slow component. 



channel is increased quadratically with p due to the in- 
creased populations of the hot electrons and hot COPs. 
Note that the scatterings among the hot electrons can- 
not account for the loss of the electronic energy 25j . At 
f ~ 1 ps, the electrons and COPs reach quasi-equilibrium 
(T e i = Tqop) after a sufficient number of scatterings, and 
the elcctron-COP composite dissipates its energy to the 
heat bath such as acoustic phonons [26|, and hence the 
p-independent energy dissipation at t > 1 ps. 

Further analysis of the spectra also supports that 
the recovery dynamics at t > 1 ps is characterized by 
T e i = Tcop an< i decay rates independent of p. First, 
we derive the decay-rate spectrum for the slower com- 
ponent 1/t2(w), which is obtained by fitting I(u>,t) at 
each energy with a double exponential function (Supple- 
mentary Information). As shown in Fig. [3th), l/r2(w) is 
independent of p and is quasi-linear to u> at u < 0.3 eV. 
Second, we quantify the spectral shape by simulat- 
ing the spectrum as I(uj,t) = a(t) J G(uj — ui',(T 2 (t) + 



a 2 R )D(uj')f(uj',T el (t))dui'. Here, G(uj,a 2 +a 2 R ) is a Gaus- 
sian with FWHM of y/a 2 + a% f(u,T el ) is the Fermi- 
Dirac function, a(t) is a scaling factor, and D(uj) is a 
DOS. The Gaussian broadening accounts for the spec- 
tral weight accumulating from lower energies, so that a 2 
becomes a measure of the number of the COPs in the 
transient. The spectra are nicely reproduced throughout 
the transient [Fig.[4ja) and Supplementary Information], 
and the fitting parameters are summarized in Fig. SJb) 
andUJc). One can see that a 2 and a for t > 1.0 ps scales 
with T e i and does not explicitly depend on p. This in- 
dicates that the number of the COPs at i>1.0ps is a 
function of T e i, i.e., Tcop = T e i, and that the spectral 
shape at £>1.0ps is determined only by T e /, which is 
the temperature of the electron-COP composite. 

Discussion. — We experimentally find that the elec- 
tronic distribution is non-thermal at t < 0.2 ps, and also 
find evidence that the COP generation is mostly accom- 
plished within this initial temporal region. This strongly 
indicates that the mechanism of the ultrafast COP gen- 
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FIG. 4: Optical-phonon broadenings during the transient, (a) 
TrPES spectra recorded under 52-mW pump (left) and the 
simulated spectra (right), (b) Fitting parameters, (c) A plot 
of the fitting parameters a 2 and a with respect to T e i for 
t > 1.0 ps. 
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cration is beyond the TTM scheme, see Fig. [5] Since 
the duration of the non-thermal electron distribution is 
longer than that in the TTM scheme, the high-energy 
electrons, which is considered to be favorable for gener- 
ating high-energy COPs [27} , indeed have more chance to 
generate COPs before they degrade into low-energy ther- 
malized electrons. Therefore, not only the largeness of 
the electron-phonon couplings but also the long durations 
of the non-thermal electron distribution may be crucial 
for understanding the ultrafast-and-efficient generation 
of the COPs in graphite. It is interesting to note that a 
break-down of TTM was also suggested in YBa2Cu307_,s 
[2e| , in which the DOS becomes vanishingly small around 
Ep in the <i-wave superconducting state, similar to the 
case for graphite. Therefore, coupled dynamics directly 
involving the non-thermal electrons may be sought in 
materials that have vanishingly small DOS around Ep, 
for example in neutral graphene [29[, in nodal supercon- 
ductors [28(, and on the surface of topological insula- 
tors [13] • Alternatively, we do not exclude the possibility 
that the hot COPs and hot electrons are co-generated by 
the pump, which may be viewed as a counterpart of the 
breakdown of adiabatic Born-Oppcnhcimer approxima- 
tion (3l|: in as much as the electrons cannot follow the 
motion of the lattice, COPs are simultaneously generated 
when the electronic excitation takes place. Whichever 
the case may be, our study reveals that there is a unique 
mechanism of ultrafast COP generation where the con- 
cept of temperatures is broken. 

Methods. — The TrPES apparatus consists of an am- 
plified Ti: sapphire laser system delivering hui = 1.5 eV 
pulses of 170-fs duration with 250-kHz repetition and a 
hemispherical analyzer (32[, see Fig. QJa). A portion of 
the laser is converted into hv4 = 5.9 eV probing pulses 
using two non-linear crystals, /3-BaB204 (BBO), and the 
time delay t from the pump is controlled by a delay stage. 
The pump and probe pulses arc p-polarized and have 
spot diameters of ~0.8 and ~0.3 mm, respectively, at the 
sample position. The intensity of the probe pulse is min- 
imized to avoid space charge effects. Multi-photon pho- 
tocmission due to the pump pulse is not observed in the 
dataset presented herein. Ep is referenced to the Fermi 
cutoff of gold and the energy resolution or is 11 meV. The 
base pressure of the photoemission chamber is lxlCP 10 
Torr, and highly-oriented pyrolitic graphite is cleaved in 
situ. 



FIG. 5: Two-temperature model (a) and beyond (b). The top 
panels show the transient change of T e i and T cop . Hatched 
area indicates temporal region where T e i cannot be defined. 
The lower panels show electron distribution functions (/) at 
some selected moments. The major difference in the COP- 
generation mechanism between the two is whether they are 
generated by the thermal electrons (a) or by the non-thermal 
electrons (b). 
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